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Abstract. 

A spherical-symmetric gamma-ray emission from the inner region of the Galaxy (at least 
up to roughly 10° in latitude and longitude) has been recently identihed in Fermi-LAT data, and 
initially associated to dark matter particle annihilations. Guided by the evidence for a high gas 
density in the inner kpc of the Galaxy correlated with a very large Supernova rate, and hence 
with ongoing cosmic-ray acceleration, we investigate instead the possibility of addressing this 
excess in terms of ordinary cosmic-ray sources and standard steady-state diffusion. We alter the 
source term, and consistently the correlated gamma-ray emissions, in the context of a template- 
htting analysis. We focus on a region of interest (ROI) dehned as: |/| < 20°; 2° < |6| < 20°, 

with I and b the Galactic longitude and latitude coordinates. We analyze in detail the overall 
goodness of the £t of our framework, and perform a detailed direct comparison against data 
examining profiles in different directions. Remarkably, the test statistic of the fit related to our 
scenario turns out to be as good as the Dark Matter one in the ROI here considered. 
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1 Introduction 

There is a remarkable agreement between models for the diffuse y-ray emission in the Galaxy 
and data from the all-sky survey from the Fermi Large Area Telescope. A fairly good match [1] 
is obtained in most region of the sky implementing only minor readjustments to a standard 
recipe based on: i) supernova remnants (SNRs) as cosmic-ray (CR) sources; ii) the steady-state 
propagation of CRs in the Galaxy as tuned on local GR measurements; and in) y-ray emitting 
targets, namely the gas and the interstellar radiation held (ISRF), as indirectly derived from 
observations at other wavelengths. 

Notable exceptions are: A systematic underestimate of the hux above the GeV in the inner 
Galactic plane [1, 2]; A hard-spectrum emission extending about 40 degrees in longitude and 50 
degrees above and below the Galactic center (GG) commonly dubbed Fermi bubbles [3] - most 
likely associated to the the signal with similar morphology previously discovered at microwave 
frequencies and dubbed WMAP haze [4]; A discrepancy between the predicted GR gradient 
and the one inferred from y-ray emissivity along the Galactic plane [5]; A roughly spherical 
excess in the inner few degrees from the GG and peaking at an energy of few GeV [6]. 

The GG region, namely the inmost few hundred pc, is one of the most challenging from the 
point of view the theoretical modelling for the diffuse y-ray emissivity. All the three ingredients 
we listed above for the standard recipe show problematic aspects. First of all, catalogues used 
for GR source distributions, mainly the observed distributions of SNRs [7] and pulsars [8, 9], 
are not optimal for the inner Galaxy. The htting functions extrapolated from such catalogues, 
and used for GR transport studies and to derive y-ray emissivities in numerical propagation 
codes such as Galprop [10-13] and DRAGON [14-16], simply go to zero at the GG and fail to 
account for the very active star forming region emerging from multi-wavelength surveys of the 
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central Galaxy. From the point of view of CR propagation, there are several features making 
it likely that transport properties in the GC region differ signihcantly from average properties 
in the Galaxy as htted to local GR data. E.g., the GG region seems to host strong magnetic 
helds, at the level of 50-100 /xG [17], much above the estimated average value within the disc or 
the GR diffusive halo, probably impacting on the structure and effective values for the diffusion 
tensor. Also, very strong convective effects have been claimed in this region [18, 19], in analogy 
to outflows from external starburst systems, as well as from the fact that at the GG there is 
a violation of the correlation between far-infrared and radio continuum luminosities observed 
for systems with in-situ energy losses of GR electrons. Finally, the density distribution and 
dynamics of molecular gas at the GG is still actively investigated, see, e.g., [20]. 

The evidence for a GG excess has been advocated based on the template-htting method. 
Extensively used for the analysis and subtraction of GMBR foregrounds, this technique has 
led hrst to the identihcation of the WMAP haze and then successfully applied to search for 
its (alleged) y-ray counterpart, the Fermi bubbles. The rationale of this procedure is to test 
whether, on top of small-scale discrepancies, the model for diffuse y-ray emission intercepts the 
correct intermediate- to large-scale morphologies, or data favour extra contributions with dif¬ 
ferent angular imprint in the region under investigation. Morphological templates are assigned 
- based on theoretical modelling and/or observations - to components connected to different 
physical emission mechanisms, such as: The component taken as the sum of the term due to 
the production and decay of vr^ plus bremsstrahlung emission, both related to the gas distri¬ 
bution in the Galaxy; the Inverse Gompton (IG) term, as correlated to the ISRF model; and 
the isotropic template for the extra-Galactic background component. Each template is then 
let free to fluctuate, independently for each energy bin, to hnd the conhguration providing the 
best £t to the data; in this way one accounts for possible spectral distortions and normalisation 
uncertainties of the different components within the theoretical model. Indications for an extra 
component can be claimed if a signihcant improvement of the overall £t is found when repeating 
the same exercise with an additional physically motivated template. 

In case of the GG the application of the template htting procedure has led to the identifi¬ 
cation of residuals with an approximately spherical morphology; such excess can be mitigated, 
and at least partially reabsorbed, introducing in the analysis an extra component which has in¬ 
deed a strong physical motivation: annihilating dark matter (DM) particles. In fact, the Milky 
Way is believed to be supported by an extended halo of DM particles, whose density should 
be maximal in the GG region; if such particles can annihilate in pairs, like in many models for 
particle-physics motivated DM candidates, they would generate a y-ray flux peaking at the GG. 
The template morphology, as well as the spectral shape found as a result of the template htting 
procedure, have been claimed to be compatible with such a DM signal, pointing to the very 
appealing scenario of Weakly Interacting Massive Particles (WIMPs) as early Universe thermal 
relics [6, 21]. Several theoretical embeddings of this picture from the DM model-building side 
were proposed [22], as well as its consistency with results from other WIMP searches has been 
tested [23]. 

On the other hand, a few competing scenarios accounting for the extra component have 
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been suggested as well. Among the proposed explanations there are: the presence a population 
of unresolved millisecond pulsars [24, 25], a bursting star-forming activity in the past near the 
GC [26], or nonthermal bremsstrahlung produced by a population of electrons interacting with 
neutral gas in molecular clouds [27, 28]. At the moment is still unclear whether there is a 
preferred scenario and further investigations are ongoing. E.g., concerning the MSP explana¬ 
tion, while there seems to be compatibility with the average spectrum of observed MSPs, the 
consistency with the number of pulsars observed by Fermi-LAT is still debated. Interestingly, 
two recent analyses aiming at an improved targeting of point sources [29, 30] suggest that the 
GC excess could be produced by a population of unresolved point sources. Also very recently, 
the burst scenario was considered in great detail in [31]; the conclusion of this work is that a 
combination of two burst events may provide a good £t of the excess spectrum, although with 
quite unnatural values of the involved parameters (e.g. extremely hard spectral indices) and 
without reproducing the observed morphology of the excess emission from the innermost degree 
or two around the GC. 

In this analysis we concentrate instead on the choice of templates for the CR-induced 
components and on the impact of such choice on the evidence for the excess. As discussed above, 
none of the ingredients of the “standard” theoretical model setting those templates are actually 
optimised for the GC environment. As a preliminary step we tested a few hypotheses on whether 
signihcant changes in the CR transport equation in GC region, including a sampling of radial 
dependences for the diffusion coefficient and anisotropic diffusion, would mimic the required 
morphological features; at this level, in agreement with and extension of results in [32, 33], we 
did not hnd suitable configurations. 

We have turned then to considering an adjustment in the choice of the CR source function, 
introducing, as a new ingredient, a steady-state source located at the GC and with a narrow 
spatial extension - few hundred pc. As already mentioned above, there are robust astrophysical 
observations supporting such new ingredient: The GC seems to harbour significant star for¬ 
mation and a large rate of Supernova explosions compared to the average value in the Galaxy; 
according to [34], the star formation rate (SFR) in the inner few degrees away from the GC is 
of order 1% of the SFR in the Galaxy, making it roughly a factor of 250 higher than the mean 
rate in the Galaxy. This should be the consequence of the following two facts: The presence 
a large reservoir of molecular gas hlling the inner part of the Galactic bulge (dubbed Central 
Molecular Zone [35] and extending up to ~ 200 300 pc away from the GC); the very peculiar 

physical properties of this environment, where the insterstellar medium appears signihcantly 
hotter, more turbulent, and more magnetized. Solid evidence for the high SFR level comes 
from infrared observations - performed with the Hubble Space Telescope - of some extremely 
dense stellar clusters in the inner 50 pc (the Central, Arches, and Quintuplet clusters). These 
structures are rich of young, very hot stars that are many times larger and more massive than 
the Sun, see, e.g., [36, 37] and references therein. Moreover, many isolated Wolf-Rayet Stars 
and O Supergiants were observed in the inner 100 pc [38]. The Supernova explosions connected 
with this relevant star formation activity are expected to accelerate a large amount of CRs. 
While this contribution has no impact on the local CR fluxes (see, e.g., the discussion on the 
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locality of CR spectra in [39]), it drives major consequences on 7 -ray emissivities in the GC 
region, as we will discuss in the rest of the paper. 

At the level of proper statistical assessments of the goodness of the £t to the data, a 
comparison between alternative theoretical models is a very hard task due to small-scale dis¬ 
crepancies which have a significant impact on statistical estimators and prevent a reliable 
comparison with the data. In order to alleviate these problems, we have adopted therefore the 
template-htting strategy. We take as reference case the £t of the data for “standard-lore” CR 
templates, rederive the improvements in the £t obtained when including the DM template, and 
then compare such rehnements to those we can obtain in the case when the new CR templates 
are adopted (different methods on how to weight £t improvements will be discussed). From 
such a comparison we will show that the two scenarios are in most respects equivalent and hence 
that the “excess” can be reabsorbed within a framework involving only a minimal picture of 
steady-state propagation of CRs. Note that the approach we follow is the appropriate one to 
realistically reproduce the energy dependence of the diffuse emission template: Indeed, while 
the Dark Matter morphology is fixed at all energies, the IC, vto and bremsstrahlung ones change 
in energy as a result of CR transport. It is therefore extremely important to incorporate the 
contribution related to the new source directly into the original diffuse emission template. 

2 Description of the method 

2.1 The template-fitting procedure 

In this section we describe the method that we use to compare our proposed astrophysical 
scenario with the DM interpretation of the GC excess. 

The starting point of the analysis is a physical model for the CR distribution in the Galaxy 
obtained with DRAGON [15, 16]. This numerical code is designed to simulate all processes related 
to Galactic CR transport (in particular: diffusion, reacceleration, convection, energy losses, 
spallation) for all CR species, from heavy nuclei to protons, antiprotons, and leptons. DRAGON 
can work in both 2D and 3D mode, and in both cases it is possible to implement anisotropic 
diffusion. The code includes the nuclear cross section database and interstellar radiation held 
model taken from the latest public version of Galprop [10]. 

We then perform a line-of-sight integral and obtain the 7 -ray skymaps from ~ 0.3 GeV to 
~ 300 GeV using GammaSky [40]. This package computes the 7 -ray emission due to: i) decay 
of neutral pions produced by collision of CRs with the interstellar gas; ii) IC scattering of CR 
leptons on the diffuse radiation held; in) Bremsstrahlung emission due to CR leptons interacting 
with the interstellar gas. Exploiting such description of the dihusion emission induced by the 
CR interactions, we model the 7 -ray sky as a linear superposition of the following templates: 

• 7 r°-decay -|- Bremsstrahlung template (produced with GammaSky). These templates are 
generated making use of the CR distributions computed with DRAGON and adopting the 
gas model taken from the public version of the Galprop CR package [10-13]. 
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• IC template (produced with GammaSky). Here we use the leptonic distributions calculated 
with DRAGON and the most recent ISRF model included in Galprop. 

• Isotropic extragalactic background template (EGB) with a spectrum taken from a recent 
analysis of the Fermi collaboration [41], 

• Point source (PS) template. This ingredient is obtained from the 4-year Point Source 
Catalog (3FGL) provided by the Fermi collaboration [42], The angular resolution of the 
Fermi-LAT instrument is taken into account smoothing the emission of each point source 
according to the point spread function (PSF) of the instrument. The PSF is modeled 
using the Fermi tool gtpsf. 


Fermi bubbles template with morphology and spectrum as derived by the Fermi collabo¬ 
ration in [43]. 

A DM template obtained integrating along the line of sight the square of a generalized 
NFW prohle: 


Pir) = Poi — 
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1 H- 
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where r corresponds is the galactocentric distance and the scale radius is fixed to 
Tg = 20 kpc. The normalization po is chosen in order to obtain a DM density at the Sun 
position p^Vq) = 0.3 GeV cm“^, with r© = 8.5 kpc. For the inner slope we take 7 = 1.26, 
as in [21]. We checked that different values of 7 do not alter the conclusion of our analysis. 
In general, we found that prohles with 7 ~ 1 — 1.6 are preferred by the data. 


We consider 5 years of Fermi-LAT data, within the event class ULTRACLEAN. We use 
the Fermi tools v9r32p5 to perform analysis cuts and to compute the exposure of the in¬ 
strument. Events with zenith angles larger than 100° are rejected. We also apply the fol¬ 
lowing selection cuts: (DATA_QUAL==1) && (LAT_C0NFIG==1) && (ABS(R0CK_ANGLE)<52) && 
(IN_SAA! =T) . The exposure maps are computed using the response function P7REP_ULTRACLEAN_V15. 
The data are organized in 30 energy bins from 300 MeV to 300 GeV, equally spaced on a log¬ 
arithmic scale. 

We work with the HEALPix tessellation scheme for all our maps [44]. We choose a resolution 
Inside = 256, corresponding to a pixel size of ~ 0.23°. Finally, the templates and the Fermi-LAT 
data are smoothed to a common angular size using a gaussian kernel. In this way all the maps 
have the same angular resolution. For this purpose we use the HEALPix routine smoothing. 

We follow the method in [45] using ftarget = 3° for < 0.6 GeV and /target = 1° for > 0.6. 

We checked that our results are insensitive to such smoothing choices, with the exception of 
the low-energy portion of the spectrum where the recipe adopted for the smoothing of the 
templates has a non negligible impact on the result, see below for further details. 

At the end of this process, we perform the template-htting analysis. We focus on a region 
of interest (ROI) defined as: 

|/|< 20 °; 2 °<| 6 |< 20 °, 


5 


with I and b the longitude and latitude coordinates. For each energy bin, we construct then 
the following Poisson likelihood function: 

-21n(£) = 2^(ei-Oi In(ei)) + 2^1n(oi!) + (2-1) 

i i 

where the index i labels each pixel inside the ROI, and we introduce the following quantities: 

• Oi is the number of observed events (the Fermi-LAT data); 

• Cj is the number of expected events. This term is written as a linear combination of the 
events corresponding to the templates described above 


e.= (2.2) 

k 

where k runs over the different templates. The counts are obtained multiplying the 
differential flux per unit of solid angle by the Fermi-LAT exposure map, the energy bin 
width and the solid angle per pixel. 


• The Fermi bubbles and the isotropic emissions are better determined with observations 
outside our ROI. For this reason, we introduce a penalty factor for their coefficients: 


Text 


E 


■A?, (1 - 


(2.3) 


where h now only refers to the Fermi bubbles and isotropic emission templates. and Ch 
are the nominal fluxes and uncertainties taken from the Fermi analysis mentioned above. 


The minimization of the function 2.1 provides the best-fit values of the coefficients in each 
energy bin. 

Concerning the treatment of point-source emission, we checked the possibility of introduc¬ 
ing a masking procedure to mitigate their role in our analysis and/or account for reconstruction 
faults. Masks of two types have considered: A sharp cut around each point source of all pixel 
within a 68% containment radius, depending on energy and sharply increasing at E <1 GeV, 
as obtained from the PSF of the instrument Fermi tool gtpsf. We also considered the soft 
mask introduced in [33], according to which the contribution of each pixel to the likelihood is 
reweigthed depending on the flux from the point sources expected in that pixel, lowering the 
impact on the likelihood of pixels whose counts are dominated by point sources. We found, 
however, that our results are not severely affected by adopting any of the two masks and hence 
disregard the masking procedure for simplicity. 
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2.2 A novel reference model 


Having set the framework, we introduce a benchmark case to examine whether, within the 
template-htting analysis, a properly modihed CR emissivity model can be considered equivalent 
to introducing the DM template. 

We choose as starting point for the CR propagation framework one of the reference cases 
in [33], labelled “Model A” in that analysis. 

This model has a standard scaling of the diffusion coefficient with rigidity (a Kolmogorov- 
type 6 = 0.33), and reasonable values for the normalization of the diffusion coefficient (5 ■ 
10^®cm^/s at 4 GV), the height of the diffusion halo (4 kpc), and the Alfven speed (32.7 
km/s); moreover, the model features a moderate convective wind gradient (50 km/s/kpc). The 
injection slope is a broken power law. For the hadrons the indexes are 1.89 (2.47) below (above) 
the reference rigidity of 11.3 GV, while for the electrons 1.6 (2.43) below (above) 2.18 GV. 

On the other hand, the model has a sharply enhanced IG component obtained by rescaling 
the normalisation of the GR electrons by a factor ~ 2 with respect to local measurements, as well 
as a ~ 40% increase in the ISRF; the synchrotron energy losses are also enhanced, adjusting the 
overall normalisation of a conventional large-scale magnetic field model for the Galaxy in such 
way that, after the rescaling, it matches the claimed value of about 50 /iG at the GG. While 
the latter assumptions make the model fail to be compatible with local GR measurements, 
we remind that the goal of the template-fitting analysis is not to find a model valid in the 
whole Galaxy, but rather to describe in detail the morphology of the GG region (reproducing 
also the local measurements, and possibily additional observables, would require instead a 
completely different method of analysis, clearly well beyond our scope; such an improvement 
will be considered in a forthcoming work). 

As also discussed in [33], the specific choice of propagation parameters is not in any way 
crucial for the template-fitting analysis. The point is that the free fioating of the individual 
templates in each energy bin always tends to readjust the input model towards a preferred 
output, carrying in the end common features, such as an enhancement in the IG component. 
Model A seems to be preferred in [33] since implies the smallest changes in normalisation 
and spectral features between input and output model (namely the minimum of the likelihood 
function is reached for the smallest departure from 1 of the coefficients introduced in Eq. 2.2). 
It is also one of the cases in which the template fitting gives a very clean indication, over a 
large span of energy bins, for a sharp improvement of the data fit when including the DM 
template; we take the challenge to focus on this sort of “worst-case scenario” to test whether 
our alternative picture can perform equally well. On the other hand, we stress again that, as 
tested on a few sample cases, changing the propagation setup would only impact on minor 
details of the fit, leaving the overall picture unchanged. 

Model A adopts as GR source function a smooth interpolation from a SNR catalogue [7] . 
As already mentioned in the Introduction, this term clearly does not include a satisfactory 
description of the Galactic bulge region. We consider then an extra steady-state source term. 
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which we model as a Gaussian (hereafter “spike”): 

Qspike Qo^Xp ^ ^2 ^ • (^•^) 

In the following, we will express the normalisation of this source in terms of TV, i.e. the 
ratio (in percentage) between the volume integral of the spike and the volume integral of 
the conventional source term. This provides a simple way to compare against observational 
constraints on the CR injection term in the center, correlated to the Star formation rate, or 
the Supernova explosion rate. 

Therefore, the most relevant parameters describing the new ingredient under consideration 
are then the spatial extension a and the normalisation M of the spike. 


IC profiles (arbitrary units) IC profiles (arbitrary units) 




Figure 1: Left panel. The IC longitudinal profiles integrated in the range 2° < |6| < 5° for 
two reference values of and for spikes with different widths (a = 200 -1- 400 pc^). We remind 
that the normalization of the spike is set in such a way to absorb most of the GC excess. We 
checked a posteriori that the chosen value is compatible with astrophysical estimates regarding 
the star formation rate in the galactic center region. Right panel. The IC profiles at 0.5 GeV 
for the spike terms with different values of a compared to the DM template. 


In our ROI the spike mainly affects the IC template, since the region containing most of 
the 7r° and bremsstrahlung emission is masked out ^ The additional contribution due to the 
spike impacts on the IC template as shown in the left panel of fig. 1 where we plot the total 

^We remark that the reference model here considered is not tuned for the Galactic plane, within 2° in 
latitude. In this region, a much more detailed treatment of the astrophysical ingredients (concerning both the 
emitting targets and the CR transport properties) would be in order; however, it is beyond the scope of our 
work. 
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IC longitudinal profile for different spike widths. On the right panel, we show instead the 
contribution to the IC emission from the spike only and compare it to the DM template. The 
plots are produced for sample 7 -ray energy values. 

It is important to notice that, while on the morphological side the two scenarios are similar, 
in the spike case we are not treating the extra ingredient as an independent term, rather we 
are correlating its spectrum to the overall IC emission. 

3 Results 



Figure 2: Spectrum of the various contributions to the totally-ray flux, pre- and post-template- 
fitting, compared to the Fermi-LAT data in the ROI characterized by |/| < 20°, and 2° < | 6 | < 
20°. The violet band superimposed to the data represents the systematic uncertainty. Left 
panel: Model A-hDM. Right panel: Model A-hspike-hDM. For Inverse Compton (light blue), 
ir^ -hBremsstrahlung (red), isotropic background emission (dark red) and Fermi bubbles (green) 
dot-dashed lines show the nominal spectrum (pre-fitting) while points and dashed lines are the 
post-fitting values. Uncertainties bands and central values for the isotropic background template 
(dark red) and Fermi bubbles (green) are taken respectively from [fl] and [43]. Magenta 
diamonds: DM contribution. The point-source template (orange dashed line) is not touched 
by the fit. The error bars on the templates are obtained from the fitting procedure. The black 
dashed (dot-dashed) line is the post (pre) template-fitting total spectrum. 

As anticipated before, let us start from the simplest realization of our idea. 

We choose a = 300 pc as a sample value for the spike width, in agreement with astro- 
physical inputs. The choice of M is based on two objective requirements:^ The value cannot be 
too large, since the SFR in the GC cannot exceed few percent of the total rate in the Galaxy; 

^We refer the reader to section 4 for a wider discussion about the impact of different values of a and N. 
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Figure 3: Residual counts obtained for the Model A (left panel, without the inclusion of DM tem¬ 
plate), for the Model A -h DM (central panel), and for the reference model described in section 3 (right 
panel, without the inclusion of DM template). See text in section 3.1 for a detailed discussion. 


moreover, we verify a posteriori that the spike emission absorbs the majority of the GC excess 
(to put it another way, it means that if we include in the fit of the reference model the DM 
template as well, the latter should be used only marginally). 

Our analysis is based on the comparison between the reference model dehned before, 
and the case in which we add to Model A a DM template, modeled according to a modified 
NFW distribution with 7 = 1.26 (see Sect. 2 for details). The inclusion of a DM template 
in the htting procedure in addition to standard astrophysical background models - like the 
Model A considered in this analysis - provided the most striking evidence supporting the 
DM interpretation of the GC excess. In [21], the addition of the DM template dramatically 
improved the fit up to an overwhelming high level of statistical significance. For this reason, it 
is mandatory to compare the performances of the spike model with those of the DM template, 
and we organize our analysis in three subsequent steps. First, in section 3.1, we present the 
result of the template £t in terms of energy spectra. Second, in section 3.2, we provide a more 
quantitative discussion based on the likelihood analysis. Finally, in section 3.3, we offer a more 
direct comparison with data based on the analysis of 7 -ray prohles along three complementary 
directions. 

3.1 Energy spectrum 

As far as the energy spectra are concerned, in fig. 2 we present the results of our template-fitting 
analysis. 

In the left panel, we show the best £t obtained for the Model A including the DM template. 
As described in the caption of fig. 2 , we show for all components both pre- and post-template 
htting values. As expected, the fit heavily uses the DM template, and the corresponding energy 
spectrum (magenta diamonds in £g. 2 ) clearly shows the familiar bump-shaped form peaked 
at = 1 - 3 GeV, in good agreement with that predicted from a WIMP with ~ 45 GeV mass 
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Figure 4; Top left panel. We compare the test-statistic (TS = — 2Alog£, we show the square- 
root of TS) of the models we consider; a positive difference between two models means that the second 
model performs better. Yellow filled circles: — 2 logTModeiA + 2 log£ModeiA+DM- Green filled circles: 
-21og£ModeiA + 21og£spike. Green empty circles: -2 log£spike+DM + 2logTspike- Top right panel. 
We compare the of the longitude profiles for the same models. Filled circles: xliodeiA “ i^ipiket 
empty circles: xliodeiA “ iTModeiA+DM- Bottom panels. The same as the top right panel, for latitude 
and radial profiles. 


annihilating, among other possibilities, into bb. For illustrative purposes, it is instructive to 
look at the residual map - i.e. the difference in counts between data and model evaluated at 
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Figure 5: Control regions used for the computation of the j-ray profiles in section 3.3. See text for 
details. 


the best-fit point - obtained without the inclusion of the DM template. We show this residual 
map in the left panel of £g. 3 for the energy window = 1 -10 GeV, where the presence of an 
excess distributed around the GC stands out in full glory. This excess is fully absorbed by the 
DM template in the £t, as shown in the central panel of fig. 3 where we present the residuals 
after template htting including the DM template. 

Let us now turn the attention to our reference model. 

In the right panel of £g. 3 we show the pattern of residuals replacing the DM template 
with the spike (A/* = 2.2%). The resulting pattern is analogous to the DM case. 

In the right panel of £g. 2, instead, we consider the possibility of having at the same 
time the spike and the DM component. The contribution of the DM template is signihcantly 
reduced, is consistent with zero (within the error bars obtained from the templated htting) in 
a large energy range, and - most importantly - gives rise to a featureless spectrum. 

In the next section we will see that precisely in this energy range, our model (without DM) 
gives a description of the y-ray emission quantitatively similar to that of the Model A-I-DM. In 
other words, most of the GG excess is absorbed by the spike, and the presence of an extra DM 
template only gives a slight improvement of the £t. 

3.2 Test statistic 

In order to scrutiny the performance of our model, let us now consider the upper left panel 
of £g. 4 where we compare, for each energy bin, the Test Statistic (TS) of the models we are 
studying. The TS is defined by TS = —2Alog£, and is computed adopting the Likelihood £ 
dehned in Eq. 2.1. 

First, we show the improvement in the TS obtained adding the DM template to the Model 
A. More precisely, the yellow dashed line with open circles in £g. 4 represents the square root of 
the difference — 2 log ^ModeiA + 2 log^ModeiA+OM (hereafter, ATS). In terms of energy spectra, 
the case Model A -|- DM corresponds to the one shown in the left panel of kg. 2. The plot 
clearly shows the improvement in the £t due to the presence of the DM template. If taken 
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comparison: latitude profiles 


Spike (III < 5°) Spike (III < 5°) Spike (III < 5°) 



-20 -10 0 10 20 -20 -10 0 10 20 -20 -10 0 10 20 
Galactic Latitude [deg] Galactle Latitude [deg] Galactic Latitude [deg] 


Model A + DM (III < 5°) Model A + DM (III < 5°) Model A + DM (III < 5°) 



-20 -10 0 10 20 -20 -10 0 10 20 -20 -10 0 10 20 
Galoctic Lolitude [deg] Golactic Lotltude [deg] Goloctic Lotitude [deg] 


Figure 6; Top row. The latitudinal profiles for our best case at several energies. From left to 
right, 400 MeV, 2 GeV, 10 GeV. Red circles (triangles): +Bremsstrahlung contribution, pre- and 
post-fitting. Blue circles (triangles): IG contribution, pre- and post-fitting. Magenta triangles: total 
model, to be compared with Fermi-LATdata. Bottom row. The same plots for Model A-hDM (DM 
contribution in forest green triangles). 


at face value, it corresponds to a statistical preference of about 15 a at the position of the 
peak. Of course we remark that this value should be taken cum grano salis, since - in addition 
to the extremely small statistical errors - data are plagued by unavoidable systematic errors 
not included in the likelihood £t. However, it is indisputable that including the DM template 
greatly improves the fit; therefore, it is crucial to compare this result with the performance of 
our spike. 

The green solid line with hlled circles represents the improvement in the £t (quantihed by 
ATS) obtained considering our reference spike model with respect to the Model A. The plot 
highlights that our scenario - without any DM contribution - performs better than the starting 
Model A, and gives a result comparable (even at the level of statistical preference) with the 


13 

























































comparison: longitude profiles 
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Figure 7: The same as fig. 6 for the longitudinal profiles. 


DM case. For completeness, we also show the TS for the combination of the spike and the DM 
template. The green solid line with open circles represents the corresponding ATS, so that 
negative values indicate a statistical preference for the addition of the DM template. In terms 
of energy spectra, this situation corresponds to the right panel of £g. 2. The TS plot shows 
that the addition of the DM template slightly improves the £t in the energy window = 1- 10 
GeV; however - as already noticed discussing the energy spectrum in fig. 2 and the residual 
map in the right panel of £g. 3 - this improvement cannot be interpreted as the evidence of a 
DM contribution since the majority of the excess was absorbed by the presence of the spike. 

Let us now pause a moment to summarize what we found. 

Looking at spectra and ATS, we clarihed that the presence of the spike - even in the simple 
realization discussed here - provides a viable astrophysical alternative to the DM interpretation 
of the GC excess. 

Nevertheless, both the energy spectra in £g. 2 and the TS in the upper left panel of £g. 4 
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comparison: radial profiles 


Spike Spike Spike 
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Figure 8; 


The same as fig. 6 for the radial profiles. 


miss a transparent comparison with the 7 -ray data thronghont the ROI. For this reason, we 
decide to get more deeply into details by analyzing latitnde, longitnde and radial prohles. 

3.3 Longitude, latitude, and radial profiles 

In order to compare models and data, we analyze the 7 -ray prohles investigating three com¬ 
plementary directions. In hg. 5 we show the control regions along the Galactic latitnde (left 
panel, red), longitnde (central panel, bine) and the projected radial distance \/P + (right 
panel, green) nsed in onr analysis. We compnte the latitnde (longitnde) prohles averaging 
on |/| ^ 5° ( 2 ° ^ | 6 | ^ 5°), and we nse 1° bins. The Galactic disk is always masked for 
\b\ ^ 2°. We show onr resnlts for the longitndinal direction in the npper-right panel of hg. 4. 
Latitnde and radial prohles are shown, respectively, in the lower-left and lower-right panel of 
hg. 4. We will come back to comment these plots in a second; to start, let us clarify some 
technical details. Our procedure goes as follows: First, we compute ~ for each direction, and 
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for each energy bin - the observed prohles of 7 -ray data. Second, using these data and the 
corresponding statistical errors, we compute the value of the function for each one of the 
analyzed models, using the best-£t coefficients obtained from the template htting procedure 
as explained in the previous section. We repeat the analysis for the Model A, the Model A -|- 
DM and our reference model. In figs. 6-8 we show, for three representative energy bins, the 
latitude, longitude and radial prohles for the reference model (top row) and the Model A -|- DM 
(bottom row). In the three aforementioned panels of £g 4, on the contrary, we superimpose 

the differences AxIm = XModeiA “ XModeiA+OM (empty circles) and Axlp^^^ = XModeiA “ xlpike 
(hlled circles). For each line, positive values indicate an improvement in the £t - compared to 
Model A - obtained including, respectively, the DM template and the spike. In each panel, 
the relative position of the two lines, versus Aygpjj^g, indicate which one of the two in¬ 

gredients, DM or spike, exhibits a better Ft. Let us now get to the heart of the matter. On a 
general ground, we will explore the three prohles separating two diherent energy ranges. On 
the one hand, the energy region = 1- 10 GeV (mid-energy region hereafter), where the DM 
template gives the most pronounced contribution; on the other one, the low-energy tail of the 
spectrum, = 0.3-1 GeV. 

o Latitude profiles. These prohles show a remarkable connection with the morphology of 
the excess. The GG excess extends far away from the GG itself, out to at least 10° in 
latitude [21, 45]. This is clear from the latitude profile shown in the central and right 
panels of the bottom row in hg. 6 where, at the reference energies E.y = 2.12 GeV and 
Ej = 10.64 GeV, the DM template gives a sizable contribution to the total count density 
up to | 6 | ~ 10°. From the central and right panels of the top row in the same hgure, we 
see that the spike gives a comparably good ht, since the role of the DM contribution is 
entirely played by the IG emission of the spike. This evidence is conhrmed in the lower-left 
panel of hg. 4 where it is clear that in the mid-energy region the spike provides a ht of the 
latitude prohles remarkably closed to one obtained by using the DM template. As far as 
the low-energy region is concerned, in the lower-left panel of hg. 4 we observe that, in the 
hrst two energy bins, the DM template produces a slightly better ht if compared with the 
spike. Let us now try to interpret this discrepancy. First of all, as already pointed out 
previously, we remind that the purpose of these plots is to oher a qualitative description of 
the comparison between models and data rather than derive strong statistical statements. 
Bearing this point in mind, we turn the attention on the left panels in the bottom and 
top rows in hg. 6 , where we zoom the latitude prohles in the second energy bin with 
central energy E.y = 0.424 GeV. The origin of the discrepancy becomes clear: while the 
DM contribution is suppressed at low energy (see also hg. 2, left panel), the steady-state 
emission of the spike, being correlated to the overall IG component, cannot show this low- 
energy dependence (at least in the simplihed model studied in this paper), and therefore 
leads to a slight overshooting in the region 3 < | 6 | < 5. 

Having said that, it is undeniable that this apparently large diherence (Xspike~XModei a+dm ~ 
50) is heavily ahected by the fact that we included in the computation of the y^s only sta- 
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tistical errors. Indeed, the systematic errors are estimated for the integrated spectra only 
(~ 5% at = 562 MeV [1]); a better nnderstanding of this uncertainty will presumably 
change the rules of the comparison, washing out, at least partially, the discrepancy. 

We therefore conclude that the spike provides a £t of the latitude prohles comparable in 
quality with the one obtained using the DM template. 

o Longitude profiles. Longitude prohles are the most delicate, since they may suffer from 
spurious emission leaked from the masked Galactic disk. As before, we start our discussion 
from the mid-energy region, and we focus our attention on the central and right panels 
of the bottom row in £g. 7. In these energy bins the bump-shaped form of the GC excess 
distinctly stands out. Even in this case, we notice that the spike provides an excellent £t 
to data, comparable in quality to the DM template. This is conhrmed by the upper-left 
panel of £g. 4, where we show the comparison between Axdm A^spike- signihcant 
discrepancy is reported, and in the mid-energy region the spike perfectly substitutes the 
DM template. 

In the low-energy region the spike improves the £t if compared with the Model A but, as 
clear from the upper-left panel of kg. 4, we find in the hrst three bins a preference (up to 
Xspike ~ XModeiA+DM ~ ^00) ^M template. In the left panels of the bottom and 

top rows in £g. 7 we show the longitude prohles in the second energy bin with central 
energy = 0.424 GeV. Glearly, the presence of the spike seems to be disfavored since 
it produces a visible overshoot of low-energy longitude prohles in |/| < 4°. As already 
stated discussing the latitude prohles, the estimate of systematic error is important to 
assess the signihcance of this discrepancy. Without them, no conclusive statement can be 
made. Most importantly, we remark that - from the point of view of transport properties 
- the GG is far from being understood; non-standard properties of dihusion, very much 
likely in such a complex environment, may completely alter the simplihed picture adopted 
so far in the description of the GG excess. Moreover, as mentioned in sec. 2, we point 
out that our results regarding the low-energy tail of the spectrum (below 1 GeV) are 
sensitive to the smoothing algorithm. In particular, if we choose the alternative strategy 
of smoothing only the model templates applying a convolution with the Fermi-LAT PSF, 
and leaving the point source, exposure, and count maps untouched, we hnd a worsening 
of the ht for F < 1 GeV. However, due to all the caveats discussed above, that energy 
region is not crucial for our analysis since it is ahected by major uncertaintes especially 
on the modelling side. Therefore, all the main results we present here are not dependent 
on these details. 

To conclude, we hnd that the spike provides a ht of longitude prohles comparable in 
quality with the one obtained using the DM template in the mid-energy region, where 
the evidence of the GG excess was considered stronger. 

o Radial profiles. Finally, we discuss the radial prohles. From the right panel in hg. 5 we 
hrst notice that the adopted binning spans a wider region of the ROI compared to the 
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previously discussed profiles. Indeed, the comparison between Axdm A^gpi^e 
lower-right panel of £g. 4 closely resembles the TS discussed in section 3.2. The spike 
provides a very good fit, comparable in quality with the DM template both in mid- and 
low-energy regions. The only exception is represented by the very hrst energy bin, but 
the same arguments outlined for latitude prohles hold true. The radial prohles in the 
top and bottom rows in £g. 8 show in detail - considering the usual benchmark bins 
at = 0.424, 2.12,10.64 GeV - the performance of the spike with respect to the DM 
template. It is clear as the modihed IC emission, altered by the presence of the spike, 
mimics the DM contribution up to \/P + ~ 10°. 

To sum up, the analysis of latitude, longitude and radial prohles enforces what already 
found in section 3.2. We conhrm that the presence of the spike in our reference model depicts 
a viable astrophysical scenario potentially able to fully explain the GC excess. 

4 Discussion 




Figure 9: TS plot for a = 100 ^ 400 pc. Figure 10; Radial for a = 100 ^ 400 pc. 

Although we presented in detail a viable working scenario, it is worth mentioning that 
several variations are possible. 

• First of all, we show that the nice agreement with the data obtained for the reference 
case, can be extended to different values of a (see £g. 9 and £g. 10). 

In particular we verihed that, going down to a = 200 pc and up to a = 400 pc (with 
a proper rescaling of Qq in Eq. 2.4), all the results of the previous section are almost 
unchanged. Interestingly, in all cases the optimal value of M does not change much. 
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We remark that spikes more extended than those considered here are in tension with 
astrophysical observations, since - as discussed in the Introduction - the evidence for a 
signihcant star-forming activity is conhned in the inner part of the Galactic bulge, and 
the SFR at i? ~ 0.5 kpc is expected to steeply decrease. 

On the other hand, for small values of a, below 200 pc, the hts to the low- up to 
intermediate-energy data signihcantly worsen, indicating an unsatisfactory description 
of the 7 -ray emission: A narrow spike, indeed, exacerbates the problem highlighted in the 
previous section with longitude prohles. 

The origin of the problem is the shape of the IC emission: once we £x the normalization 
of the spike in order to match the 7 -ray emission at energies above ~ 1 GeV, the IG 
emission for a narrow spike is signihcantly steeper than larger values of a, as shown in 
£g. 1. Such morphology is disfavored by the low energy data. 

• In order to provide a more realistic description of the inner Galaxy, one should take into 
account complex interplays in the GG environment impacting on GR tranport properties. 

In particular, we know that the regular Galactic magnetic held can be decomposed into 
several components: besides a well-known pattern lying on the Galactic plane, in a recent 
detailed analysis [46] the authors pointed out a X-shaped component extending from the 
GG up to a Galactocentric radius R ~ 4-75 kpc and to z ^ 2-73 kpc. Since GR dihusion is 
expected to be anisotropic, and since - according to the Quasi-Linear Theory (QLT) valid 
in the low-turbulence limit - the parallel dihusion should dominate over the perpendicular 
one, we expect an efficient GR escape due to parallel dihusion in the vertical direction. 

We checked several values of the ratio between the perpendicular and parallel dihusion 
components of the dihusion tensor {Dzz/Drr) for each value of a, but we did not hnd 
any relevant and solid improvement in the quality of the ht valid in the whole energy 
range; We postpone to a more detailed forthcoming work a careful study regarding the 
impact of this ingredient. 

• Although we chose Model A as a reference case, we verihed that our results are not 
strongly dependent on the GR transport model. 

The same trends we discussed in detail in the previous section emerge if other GR prop¬ 
agation models are adopted. In particular, we considered the conventional KRA model 
described in [39], and Model F in [33]. In both cases, we hnd that the spike plays a crucial 
role when the template-htting algorithm is applied, and all the results presented in the 
previous section still hold. 

5 Conclusions 

This paper investigates the possibility of addressing the 7 -ray Galactic center excess, as pointed 

out in the literature and possibly associated to a DM annihilation component, in terms of 
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ordinary SNR sources and standard steady state CR diffusion. A new reference model is 
introduced including a population of SNRs in the central region of the Galaxy; while such 
term was previously overlooked in CR propagation modelling, it is strongly motivated based 
on the observed level of star formation rate and supernova explosions in the central region of 
the Galaxy. The main impact of this extra source in the Galactic center analysis is to enhance 
the IC emissivity, introducing a morphological imprint analogous - but different when coming 
to details - to DM pair annihilations in a spherical NFW-like density profile. 

Differently from previous proposals for alternative scenarios, we introduced the new in¬ 
gredient in the context of the template-fitting algorithm, treating the new contribution as 
correlated to the conventional IC emission. We analyzed in detail the overall goodness of the 
fit of our framework, and produced a more direct comparison against data examining profiles 
in different directions. Remarkably, the test statistic of the fit related to our scenario turns out 
to be as good as the Dark Matter one in the ROI here considered. 

Our result is solid against different point-source masking technique, standard choices of 
the CR transport properties, and holds for a central source width ranging from ~ 200 to ~ 400 
pc. 

The additional source proposed here provides extra 7 -ray emissivity also at energies larger 
than the GeV range analyzed in this work, although with a progressively smaller angular extent 
given that, due to the increase in the energy loss efficiency, electron diffusion is reduced on a 
shorter scale. This prediction can be tested in the near future, in the perspective of comparing 
low- and high-energy measurements by experiments such as ASTROGAM [47] and GAMMA- 
400 [48]. 

Although the template-fitting method is very powerful in addressing morphological agree¬ 
ments or flaws of a given theoretical model, there are some drawbacks as well in this approach: 
Since the free-floating of the templates is performed independently on each energy bin, the 7 -ray 
spectra corresponding to the different components can be seriously altered, resulting sometimes 
in little control on the way the physical properties of the GR transport model are modified. 
Motivated by these considerations, as a further development of this work, we plan to exploit an 
alternative approach in which the original features of a physical propagation model - supposed 
to provide also a good fit of the local observables - are kept under control by the data fitting 
procedure. 
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